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ABSTRACT 

We consider in detail the spectral energy distribution (SED) and multi-wavelength variability of 
NGC 5548. Comparison with the SEDs of other AGNs implies that the internal reddening of NGC 5548 
is E(B-V) = 0.17 mag. The extinction curve is consistent with the mean curve of other AGNs found 
by Gaskell & Benker, but inconsistent with an SMC-type reddening curve. Because most IR emission 
originates exterior to the broad- line region (BLR), the SED seen by the inner BLR is different from 
that seen by the outer BLR and from the earth. The most likely BLR covering factor is ~ 40% and it 
is not possible to get an overall BLR covering factor of less than 20%. This requires that the BLR is 
not spherically symmetric and that we are viewing through a hole. Line-continuum variability transfer 
functions are consistent with this geometry. The covering factor and geometry imply that near the 
equatorial plane the BLR covering approaches 100%. The spectrum seen by the outer regions of the 
BLR and by the torus is thus modified by the absorption in the inner BLR. This shielding solves the 
problem of observed BLR ionization stratification being much greater than implied by photoionization 
models. The BLR obscuration also removes the problem of the torus covering factor being greater 
than the BLR covering factor, and gives consistency with the observed fraction of obscured AGNs. 
The flux reduction at the torus also reduces the problem of AGN dust-reverberation lags giving sizes 
smaller than the dust-sublimation radii. 

Subject headings: galaxies: active — galaxies: quasars: general — X-rays:galaxies — black hole physics 
— accretion: accretion disk — dust, extinction 



1. INTRODUCTION 

Because of the difficulty of observing in the extreme 
ultra-violet (EUV), the true shape of the continuum of 
AGNs has long been a mystery. If emission lines arise 
from photoionization, then the equivalent width of the 
lines gives information on the number of ionizing photons 
(Zanstra 1931, Osterbrock 1989) and their total energy 
(Stoy 1933). Early photoionization model calculations 
of the broad- line region (BLR) spectrum (e.g., David- 
son 1972, MacAlpine 1972, Shields 1972) assumed cov- 
ering factors (Q/4tt) and continuum shapes that could 
account for the observed strengths of the emission lines. 
However, the first spectrophotometry of the rest-frame 
UV continuum in high-redshift AGNs showed an appar- 
ent turn down of the UV spectrum (Oke 1970, 1974). 
IUE satellite observations of lower-redshift AGNs con- 
firmed this (Green et al. 1980). Green et al. reported 
that the a « 0.5 spectral slope (F v cx v~ a ) from the 
optical to the UV steepened to a « 2.3 at wavelengths 
shorter than A1200. MacAlpine (1981) showed that there 
was a severe "energy budget" problem because such a 



1 Present Address: 
sity of Texas, Austin 



Department of Astronomy, Univer- 
TX 78712-0259. Electronic address: 
gaskell@astro.as.utexas.edu 

2 Present Address: Department of Astronomy, New Mexico State 
University, Las Cruces, NM 88003-0001. Electronic address: ek- 
limek@nmsu.edu 

3 Electronic address: lsn@kzn.ru 

4 Visiting scientist, Department of Physics and Astronomy, Uni- 
versity of Nebraska 



continuum needed S1/4-7T > 9 to match the observed He 
II A4686 equivalent width. The observed ionizing con- 
tinuum was thus providing an order-of-magnitude fewer 
ionizing photons than were needed. 

The problem was compounded when spectroscopy of 
the Lyman-limit region showed that rest-frame Lyman 
continuum absorption was rare (Osmer 1979). Smith et 
al. (1981) concluded from observations of the Lyman- 
limit region that Q/Att < 0.15 and further analysis by 
MacAlpine (1981) showed that Q,/An « 0.05 was most 
likely. Netzer (1985) pointed out that there is a seri- 
ous energy-budget discrepancy for other lines. The HST 
composite spectrum of Zheng et al. (1997) and the multi- 
wavelength spectra of a sample of PG quasars (Laor 
1997) confirmed the steepening of the observed contin- 
uum at shorter wavelengths. Photoionization models 
(e.g., Shields & Ferland 1993; Korista et al. 1998; Goad 
& Koratkar 1998; Kaspi & Netzer 1999; Korista & Goad 
2000) consistently required covering factors an order of 
magnitude higher than allowed by the limits on Lyman 
limit absorption. 

An important question is whether the continuum we 
see is the same as that seen by the emission lines. The 
continuum could be intrinsically anisotropic, or the con- 
tinuum we see could be modified by extinction from dust 
grains along the line of sight, or by absorption lines. 
MacAlpine (1981) suggested that a marked upturn in 
dust extinction below A1200 could be causing the appar- 
ent turndown in the continuum we see. 

Green et al. (1980), however, felt that the turndown 
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below A1200 was probably intrinsic, and that extinction 
by dust grains was probably not causing it because of the 
lack of effect on the spectrum longwards of A 1200. Mc- 
Kee & Petrosian (1974) had previously argued against 
the presence of dust in AGNs because of the lack of both 
the A2175 dust feature and the curvature in the UV spec- 
tra that dust should produce. The apparent match up 
between UV spectral indices and the X-ray region (Laor 
1997) argued against significant extinction. 

Various other solutions have been offered to the energy- 
budget problem. Collin-Souffrin (1986), Joly (1987), 
and Dumont, Collin-Souffrin, & Nazarova (1998) favored 
non-radiative heating of the BLR as the solution to the 
energy-balance problem. Binette et al. (1993) suggested 
the clouds were being heated by a diffuse continuum very 
close to them. Korista, Ferland, & Baldwin (1997) sug- 
gested the the UV-EUV SED might be double peaked. 
Maiolino ct al. (2001c), on the other hand, suggested 
that the solution to the conflict between the lack of Ly- 
man limit absorption and the high covering factors im- 
plied by photoionization models was that the BLR was 
not covering the source uniformly. 

In this paper we investigate the energy-budget question 
in the well-studied AGN NGC 5548 and investigate the 
implications for the BLR and torus structure. NGC 5548 
is an important test case for investigating the energy- 
budget issue (Rokaki et al. 1994; Dumont et al. 1998) 
since NGC 5548 has been particularly well studied at all 
accessible wavelengths, its variability behavior in most 
spectral regions is well known, and the line-continuum 
transfer functions, which give information on the BLR 
size and radial distribution, are better known than for 
any other AGN. We focus in this paper on the period 
2449115 - 2449130 when NGC 5548 was intensely mon- 
itored by the Hubble Space Telescope. We adopt a dis- 
tance, R, of 76 Mpc to NGC 5548. This is based on an 

H = 73 km/s/MpC, flmatter = 0.27, ^vacuum = 0.73 

cosmology and the local velocity field model of Mould et 
al. (2000). 

2. THE SPECTRAL ENERGY DISTRIBUTION OF 

NGC 5548 

Many different spectral energy distributions (SEDs) 
have been considered when modelling NGC 5548. A 
number of these SEDs are illustrated in Fig. 2 of Du- 
mont ct al. (1998), in Chiang & Blaes (2003), and Fig. 1 
of Steenbrugge et al. (2005). It can be seen that the SEDs 
differ by up to factors of three in the region below the 
Lyman limit and elsewhere. We therefore start off by re- 
constructing the probable SED of NGC 5548 around JD 
2449120 in detail. 

Like all AGNs, NGC 5548 varies. It was in fact the 
very first AGN for which variability was discovered (see 
historical discussion in Gaskell & Klimek 2003). This 
variability has to be allowed for in constructing the con- 
tinuum shape. As far as possible, one needs simulta- 
neous observations of the continuum at different wave- 
lengths. However, such observations are not available 
for all wavelength regions, so, where there were no si- 
multaneous observations, we have estimated the fluxes 
by scaling observations from other epochs. For deter- 
mining the appropriate scalings we have used other sets 
of quasi-simultaneous data. These scalings are necessary 
for estimating the SED at a given epoch, but they should 



TABLE 1 

Observed and Corrected Continuum Fluxes 



Region 


Log(v) 


Log(!AF„) a 


MW de-red 
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100 fim 


12.48 


-10.14 






60 fim 


12.70 


-10.29 






25 fim 


13.08 


-10.00 






12 fim 
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-10.03 
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-10.32 
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14.13 


-10.10 






1.65 fim 


14.27 


-10.03 






1.25 fim 


14.38 


-10.06 






5200 A 


14.76 


-10.24 


-10.21 


-10.13 


4820 A 


14.79 


-10.25 


-10.22 


-10.10 


4650 A 


14.81 


-10.27 


-10.23 


-10.10 


2250 A 


15.12 


-10.29 


-10.20 


-9.72 


2050 A 


15.17 


-10.29 


-10.20 


-9.69 


1820 A 


15.22 


-10.24 


-10.17 


-9.63 


1490 A 


15.30 


-10.24 


-10.16 


-9.59 


1360 A 


15.34 


-10.21 


-10.13 


-9.54 


1179 A 


15.41 


-10.23 


-10.12 


-9.53 



a vF v in ergs s 1 cm 2 



also be useful for modeling spectral variability at differ- 
ent epochs. 

From the 73-day period JD 2449061 to 2449134 
NGC 5548 was simultaneously monitored with the Hub- 
ble Space Telescope, the IUE satellite, and many optical 
observatories (Korista et al. 1995). From the period of 
HST coverage we have selected JD 2449115 - 2449130 
when NGC 5548 was in a high state. During this pe- 
riod data were available only for some of the optical and 
UV continua and emission lines, so in the following sub- 
sections we describe, in order of construction, how the 
continuum in each waveband was obtained. 

2.1. Optical and UV 

In Table 1 we give our measured optical continua from 
the International AGN Watch (I AW) optical spectra (see 
Wanders & Peterson 1996 5 ). Three relatively line- free 
sections of the spectra between major emission lines were 
taken to represent the continuum. These sections were 
defined by the observed wavelength ranges: 4640-4660 
A, 4810-4830 A, and 5190-5210 A. We measured average 
UV continuum points from the HST spectra (Korista et 
al. 1995) during the high state at AA 1170, 1360, 1490, 
1820. 2050, and 2250 A. We used observed wavelength 
ranges of 1160-180 A, 1344-1374 A, 1478-1501 A, 1804- 
1834 A, 2022-2204 A, and 2225-2241 A. 

To obtain fluxes in other passbands, when there were 
no simultaneous observations, we scaled observations to 
our optical or UV fluxes. There are many broad-band 
optical photometric observations available but these in- 
clude host-galaxy contamination. We need to correct for 
different host galaxy contributions and differing pass- 
bands. The I AW database gives monochromatic fluxes 
at 5100A standardized to a 5 x 7.6 arcsecond spectro- 
scopic aperture. Romanishin et al. (1995) discuss how 
to estimate the nuclear flux of NGC 5548 from direct 
images and they give the transformation between the 

5 Note that all the line and continuum fluxes in their paper are 
too low by a factor of two. 



AGN Broad-Line Region Geometry 



3 



5100A fluxes and V-band fluxes. They give estimates of 
AGN-free galaxy magnitudes as a function of photomet- 
ric aperture and give the host galaxy light contribution to 
the standard I AW spectroscopic aperture. They estimate 
the galaxy flux in the standard photometric aperture to 
be 2.53 x 1CP 15 ergs s _1 cm~ 2 A^ 1 , while the galaxy con- 
tribution in an 8-arcsecond radius photometric aperture 
is 5.81 x lCT 15 ergs s _1 cm~ 2 A -1 . 

Although for our modeling below we consider the host- 
galaxy-free AGN optical fluxes, for comparison with 
other wavelengths, our observed fluxes (column 2 of Ta- 
bic 1) still include the galaxy contribution since this is 
included by other observers. In section 4 wc subtract 
host galaxy contributions from all wavebands. 

2.2. Infra-red 

Even though infra-red radiation is not relevant in pro- 
viding ionizing photons, it does contribute to the heating 
of the broad line region and hence to the total energy 
budget. As we will discuss below (see section 2.3), in or- 
der to know the IR flux received by the BLR clouds one 
needs to know the distance of the IR emitting regions 
from the center. It will also be important to know the 
inner radius of the dusty torus because, as proposed by 
Netzer & Laor (1993) and Laor (2007), we will find that 
this is probably the outer edge of the BLR. 

To our knowledge, there were no simultaneous infrared 
observations made at the time of the HST observations. 
We can, however, estimate the IR flux from the relation- 
ship between the IR and optical at other times. From a 
study of 41 Seyfert galaxies, Glass (2004) reported that 
all but two showed variability, with K-band (2.2 fim) am- 
plitudes in the range < 0.1 to > 1.1 mag. The shortest 
timescale he found for detectable changes was one week. 
Glass found that, in general, the near-IR flux varies in 
the same way as the optical, UV, and X-ray flux, but 
with a time delay. For NGC 5548 he reports observa- 
tions in the J, H, K, & L bands made on 12 nights over 
the three-year period 1988 - 1990. Lyutyi & Doroshcnko 
(1993) give UBV photometry during the same period. 

The delay of the IR fluxes relative to the optical and 
UV is substantial so the IR flux depends on the optical 
flux at an earlier time. Furthermore there will be an ad- 
ditional delay in IR radiation from the dust getting back 
to the BLR clouds. The lag of the IR flux was not mea- 
sured for NGC 5548 in the high state we consider, but we 
can estimate the lag by scaling from the measurements of 
the lag by Suganuma et al. (2004) when NGC 5548 was 
in a low state and independently by scaling from other 
objects. 

Glass (2004) gives estimated delays and 3.4/im lumi- 
nosities for 17 AGNs. In Fig. 1 we show the relationship 
between the IR lag, r, and luminosity. A least squares 
fit gives 



log(r) =0.481og(L 3 .4^ 



9.14 



(1) 



This correlation is completely consistent with the ex- 
pected r oc L 5 relationship for a constant dust subli- 
mation temperature. Minezaki et al. (2005) and Sug- 
anuma et al. (2006) find a similar correlation between 
r and the V-band absolute magnitude. For NGC 5548, 
L3. 4pm = 22.81, so Fig. 1 implies a delay of <~ 84 days. 

Suganuma et al. (2004) monitored NGC 5548 in the 




Log(L 3 . 4(lm ) 

Fig. 1. — IR lags, r (in days), as a function of L-band flux (log 
WHz ) for 17 AGNs. The lags are the responsivity-weighted lags 
of the J and L band fluxes relative to the U band fluxes. The line 
is a least-square fit. Data from Glass (2004). 

IR and optical bands from JD 2451950 to JD 2452825, 
when NGC 5548 was in a low state. They measured 
time delays between the K band and the V band of 48 ± 
3d and 47 ± 6d. Even though the mean V-band flux of 
NGC 5548 decreased by a factor of about 2.8 between 
the earlier epoch we are considering, their lags can be 
used to check our estimate if we consider the change in 
the level of activity. Due to dust sublimation, an increase 
in luminosity leads to an increase in the effective torus 
radius, thereby increasing the lag between the IR and the 
optical according to the relationship r oc ^JL (Oknyanskij 
et al. 1999.) Thus, when the measured lag is scaled by the 
ratio of the average luminosities during the two epochs, 
we should get a lag consistent with our estimate. Scaling 
of the Sugnuma et al. lag gives 79 d, which is in good 
agreement with our earlier estimate of 84 d from the 
general IR radius luminosity relationship, and within the 
error bars of the Suganuma et al. lags. 

By comparing the J-band (1.25 (im) fluxes of Glass 
(2004) with the average V magnitudes 80 days earlier, 
from Lyutyi & Doroshenko (1993) we get the approxi- 
mate relationship: 



J = 0.4V + 6.2 



(2) 



To estimate the IR continuum seen by the BLR from 
the optical continuum we need to look at the optical 
continuum, not at a time r earlier, but at a time ~ 2r 
earlier. This is because the UV/optical continuum heat- 
ing the dust takes <~ r days to get to the dust torus, 
and the reprossed IR takes another ~ r days on aver- 
age to get back to the BLR. Because of the larger size 
of the IR-emitting region, we took the average optical 
flux (from Korista et al. 1995) between 2r — 50 days be- 
fore and 2r + 50 days before our high and low states. 
After appropriate scaling of the Korista et al. fluxes to 
the magnitudes of Lyutyi & Doroshenko (1993), we ob- 
tained J magnitudes of 11.63 and 11.66 for our high and 
low states. As would be expected, the IR flux does not 
follow the more rapid UV/optical changes. 

Choloniewski (1981) showed that, in the optical, the 
shape of the variable component of AGNs remains the 
same as the flux level changes (see also Winkler et al. 
1992). The analysis of Glass (2004) shows that, for any 
given AGN, the colors of the variable component of the 
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IR emission are also usually independent of the level of 
activity and, for his whole sample, he finds that the col- 
ors of the variable components fall within moderately 
narrow ranges. We have therefore used our calculated 
J magnitude and the average IR colors of Glass to get 
estimated H, K, and L magnitudes of 10.77, 10.05, and 
8.76 respectively for the low state, and 10.74, 10.02, and 
8.73 for the high state. These fluxes fall well within the 
historical ranges of fluxes given in the literature. 

For longer IR wavelengths we have taken measure- 
ments from the literature. Potential problems with 
these measurements include inhomogeneity due to vary- 
ing data quality and photometric apertures used, and 
source variability. We excluded measurements that were 
taken using apertures larger than 15" diameter, as these 
measurements contain significant additional off-nuclear 
infrared emission from the host galaxy. For the M 
(4.8 /im) waveband we averaged the fluxes of McAlary, 
McLaren, & Crabtree (1979) and McAlary et al. (1983), 
and for the N (10 /im) waveband we averaged the fluxes 
of Kleinmann & Low (1970), Rieke & Low (1972), and 
Rieke (1978). There is observational evidence that for 
these longer IR wavelengths non-simultaneity is not a 
problem. Neugebauer et al. (1989) found that while 
AGNs show variability in J, H, K, and L, there was no 
variability at 10 /im. Similarly, Edelson & Malkan (1987) 
found that there was no evidence for variability in the 
IRAS passbands for non-blazar AGNs. 

The 12, 25, 60, and 100 /an IRAS observations (Moshir 
et al. 1990) have a low spatial resolution and no variabil- 
ity information is available. They refer to a region much 
further out in the galaxy than the near-IR torus. 

2.3. IR Continuum Seen by the BLR 

A large fraction of the IR continuum in an AGN is 
clearly not coming from the central source, but from a 
torus well outside the BLR. The ratio of IR to optical/UV 
that a distant observer sees is not the same as what a 
BLR cloud sees. 

Calculating the X-ray, UV, and optical flux at the BLR 
clouds is straight forward. We assume that the region 
emitting this flux is smaller than the BLR and located 
inside it. We can then obtain the flux by simply mul- 
tiplying the flux observed at the earth by (R/r) 2 where 
R = 9.05 x 10 10 Id is the distance to NGC 5548, and r is 
the distance of the BLR clouds from the center. For the 
IR, for which a substantial part of the emission observed 
from the earth arises from emission from the torus out- 
side the BLR, the assumption of a compact source within 
the BLR is incorrect, and we will overestimate the flux 
at the BLR if we just scale the flux in the same way as 
the X-ray, UV, and optical. Instead, the IR flux at the 
BLR clouds depends on the surface brightness of the IR 
emitting regions and the solid angles they subtend. 

We roughly estimated the relative dilution of the IR 
continuum at the BLR clouds with the following scheme. 
We divided the radiation received by a BLR cloud at ra- 
dius r c into two components: the radiation originating 
from r < r c and the radiation originating from r > r c . 
We assumed that the flux of radiation originating at 
r < r c could be scaled from the flux received at the earth 
by (R/r) 2 , where R is the distance to NGC 5548. If we 
consider a uniform thin spherical emitting shell of radius 
r shell! the radiation density falls off as r~ 2 for r > r s heii 



while it is constant for r < r s heii (i.e., within the shell). 
This still remains a reasonable approximation if the emit- 
ting matter is in a ring rather than a shell. Therefore, 
for radiation originating outside r c the flux was scaled 
by (R/ r s h e u) 2 . Thus the IR radiation contribution from 
outside the BLR was reduced by (r c /R s heii) 2 relative to 
the contribution from inside the BLR. 

To estimate the fraction of IR radiation from within 
r c we assumed there was a minimum accretion-disk tem- 
perature of 5000 K. Any IR flux less than that produced 
by the spectrum of a 5000 K black body was assumed 
to come from outside r c and reduced by (r c / Rsheii) 2 . 
In a standard thin accretion disk (Pringle & Rees 1972, 
Shakura & Sunyaev 1973) the temperature, T falls off 
with radius as r~ - 75 . If there is reprocessing of radiation 
from the inner regions to the outer regions then the tem- 
perature gradient is shallower (e.g., T cx R~ ' 50 for the 
"slim disk" of Abramowicz et al. 1988). If dust is in ther- 
mal equilibrium with radiation from a compact source, 
Tdust oc R~ 50 as well, so we assume that T cx R~ 50 . 
This implies that the effective radius emitting a given 
wavelength is R e ff oc A 2 . The reduction of the IR ra- 
diation is thus cx A 4 . The contribution from outside the 
BLR can thus be ignored at all but the shortest wave- 
lengths. 

Because of this effect, BLR clouds at different radii re- 
ceive different long-wavelength SEDs. The outer BLR 
is close to the inner edge of the torus and experiences 
a K-band flux similar to that given by a simple inverse- 
square-law scaling; the inner BLR is an order of magni- 
tude closer and receives 100 times less relative K-band 
flux from the dust than would be expected by simple 
scaling. In our photoionization modeling we therefore 
use different SEDs at different radii. 

2.4. Hard X-rays 

2.4.1. Relationship Between Hard X-ray and Optical/UV 

Flux Levels 

There were no simultaneous hard X-ray observations 
during the HST campaign, so we need to estimate the 
hard X-ray flux from the optical and UV. Uttley et al. 
(2003) show that the 2-10 keV flux measured with RXTE 
is well correlated on long timescales (months) with the 
optical flux at A5100. However, on short timescales (day- 
to-day), the X-rays in NGC 5548 can vary a lot while the 
optical varies little (see Uttley et al. 2003 and Gaskcll et 
al. in preparation). We compared the RXTE 2-10 keV 
fluxes from Uttley et al. (2003) with the 5100A fluxes 
(from the IAW database) measured within ±7 d. The 
median separation time between the X-ray and optical 
observations is 0.6 d. As can be seen in Fig. 2, there is a 
direct proportionality, after subtracting the galaxy con- 
tribution as estimated by Romanishin et al. (1995). We 
found that, apart from short-term events on the timescale 
of a month or two, the average X-ray/optical ratio was 
constant in time and showed no correlation with the op- 
tical flux level. The X-ray flux, in units of 10 -11 ergs s _1 
cm~ 2 , can be predicted from the optical by the equation 

F 2 -wkeV = 0.564 x F X5W0 , (3) 

where the iVsioo is in units of 10 -15 ergs s _1 cm~ 2 A -1 . 

Clavel et al. (1992) (see their Fig. 4) found that the 
2-10 keV flux observed with Ginga, was also directly 
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Fig. 2. — The relationship between the 2-10keV flux (in units 
of 10~ n ergs s _1 cm -2 ) and the host-galaxy subtracted optical 
flux (in units of 10 -15 ergs s _1 cm~ 2 A -1 ) for near simultaneous 
observations. The straight line is the relationship i^-lOfceV = 
0.564 x F A5100 

proportional to the 1350A flux: 

F 2 -10keV = 0.114 X F A i 350 . (4) 

In our high state, F\ 5W0 = 6.49 x 10~ 15 after host- 
galaxy subtraction, and i*Ai35u = 44.8 x 10" 15 . From 
equations (3) and (4) these predict F 2 -iokeV = 3.7 x 
10~ n and 5.1 x 10 -11 ergs s _1 cm" 2 respectively. We 
adopt a mean of 4.4 x 10 -11 ergs s _1 cm~ 2 . 

2.4.2. Shape of the Hard X-ray Spectrum 

As in most AGNs, the shape of the hard X-ray spec- 
trum of NGC 5548 can be described by a power law 
(F v oc v~ a ) with a spectral index a <~ 0.7. In NGC 
5548 a depends on the flux level. Markowitz, Edel- 
son, & Vaughan (2003) show, from RXTE data, that for 
NGC 5548 the degree of hard X-ray variability is con- 
stant as a function of energy on short timescales but 
decreases for higher energies on long timescales. Chi- 
ang et al. (2000) and Papadakis et al. (2002) find a good 
correlation between a and F 2 ~iokeV- From Fig. 4 of Chi- 
ang et al. and our estimated F 2 -iokeV fluxes we get a = 
0.76. This is within the range found by Magdziarz et 
al. (1998) and also consistent with the trend found by 
Walter & Courvoisier (1990). 

The hard cut-off in the 7-ray region is well constrained 
by OSSE spectra from the Compton Gamma Ray Obser- 
vatory. Magdziarz et al. (1998) get a hard energy cutoff 
of 120 keV for 1991 and 1993 OSSE observations. 

2.5. Extreme Ultraviolet and Soft X-ray Region 

The most energetically important region of the spec- 
trum of AGNs is the extreme ultraviolet (EUV) and soft 
X-ray spectral region. This is also the region responsible 
for most of the ionization of the BLR. Unfortunately, it 
is very hard to observe AGNs in the EUV and soft X-ray 
spectral region because of heavy atomic absorption. 

Walter et al. (1994) present ROSAT soft X-ray obser- 
vations of several AGNs, including NGC 5548. For all 
objects they find a strong soft X-ray excess above the ex- 
trapolation of the higher energy X-rays to lower energies. 
They find that the soft X-ray spectral shape is very sim- 
ilar in the different sources even though they range over 



four orders of magnitude in luminosity. The soft X-ray 
excess can be described by a bump with a similar cutoff 
energy of kT ~ 200 keV for sources with a wide range 
of luminosity. For NGC 5548 itself, the effective tem- 
perature of the soft excess appears to remain constant 
at kT ~ 400 eV (see Fig. 7b of Magdziarz et al. 1998). 
More recent XMM-Newton and Beppo-SAX observa- 
tions (Pounds et al. 2003) also show a clear soft excess 
below 700 eV which, after allowance for complex absorp- 
tion, probably smoothly extends to 2 keV in agreement 
with the Walter et al. (1994) spectrum. There is debate 
over whether the derived temperature of the soft excess in 
AGNs is real, or whether the apparent exponential cut- 
off is due instead to strong relativistically smearing of 
a partially-ionized absorption (Gierlinski & Done 2004), 
but the fundamental nature of the soft excess does not 
affect the continuum the broad-line region sees, so long 
as the BLR sees the same continuum we see. 

There was no simultaneous monitoring of the soft X- 
ray spectrum during the HST/IUE campaign. We must 
therefore try to estimate the EUV/soft-X-ray flux from 
either the UV or the hard X-ray components. To do 
this we need to consider the nature of EUV variability 
and the relationship of EUV and soft X-ray variability to 
variability of the continuum at higher and lower energies. 

2.5.1. Variability of the Soft Component 

Our knowledge of variability of the EUV/soft X-ray 
flux of NGC 5548 has come primarily from monitoring 
with the ROSAT satellite, covering 0.1 - 2.4 keV, and 
the EUVE satellite, giving effectively measurements at a 
single energy of <~ 0.2 keV. Walter & Courvoisier (1990) 
found that as NGC 5548 varied, although the ratio of 
the soft excess to the higher energy power law compo- 
nent changed, the shape of the soft X-ray spectrum did 
not. EUVE count rates cannot be readily converted into 
flux units because the conversion depends on the spec- 
tral shape and hydrogen column density (see Marshall et 
al. 1995), but given the apparent constancy in the spec- 
tral shape, they are important for showing the relative 
variability at <~ 0.2 keV. 

Observations reveal two important things about the 
soft excess. The first is that the soft excess varies most. 
The ratio of changes in the hard and soft X-ray com- 
ponents during the ROSAT monitoring of Nandra et al. 
(1993) is consistent with the soft excess varying more, 
and the ROSA T monitoring of Done et al. (1995) shows 
that variability is primarily at energies of less than 0.4 
keV (see their Fig. 2). From subsequent simultaneous 
EUVE and X-ray monitoring of NGC 5548, Chiang et al. 
(2000) found that the amplitude of the EUVE variability 
was 15 - 20% greater than that of the RXTE variability. 

The second important thing is that, on short 
timescales, the soft excess leads the variability. Kaas- 
tra & Barr (1987) discovered from EXOSA T monitoring 
of NGC 5548 that the soft X-rays led the hard X-rays 
by about 1.5 hours. Walter & Courvoisier (1990) found 
a similar leading of about 1 hour. Chiang et al. (2000) 
found that the EUV (0.2 keV) led the 0.5-1 keV ASCA 
band in NGC 5548 by 0.15 ± 0.05d, and led the RXTE 
2-20 keV band by 0.45 ± 0.1 Id. These small delays are 
unimportant for our analysis, but the sign of the delay 
together with the greater amplitude of variability of the 
soft excess implies that short-term X-ray variability is 
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Fig. 3.— UV fluxes at 1350A, observed with the IUE satellite 
versus EUVE 84.A fluxes (from Marshall et al. 1997) averaged over 
±ld. The solid line is a least-squares fit. The dashed line represents 
a constant ratio of EUV to UV flux. 

driven by variability of the soft excess. The width at 
zero intensity of the EUV/hard X-ray cross correlation 
function of Chiang et al. (2000) is ±ld. The size of the 
EUV emitting region is thus < 1 lt.d. The EUVE light 
curve of Chiang et al. (2000) shows there are sometimes 
changes of up to a factor of two on a timescale of an 
hour. We can thus be confident that the EUV to soft 
X-ray emission originates well inside the BLR. 

2.5.2. Correlation of Soft Excess with Lower Energies 

Although there were no simultaneous soft X-ray spec- 
tral measurements during the 1993 HST campaign, Mar- 
shall et al. (1997) made EUVE observations which over- 
lap some of the UV observations, and despite the prob- 
lem of converting EUVE counts to fluxes, the amplitude 
of EUV/soft X-ray variability can be estimated from the 
EUVE observations since the EUVE samples part of the 
ROSAT band at <~ 0.2 keV. This monitoring is impor- 
tant for showing the relative change in the flux of the soft 
excess as a function of the UV flux. 

The EUVE observations over the period JD 2409104 - 
2409112 (shortly after the high state we are considering) 
overlap the HST monitoring, and other EUVE observa- 
tions overlap the IUE monitoring. The amplitude of the 
EUV variability is about <~ 5 times greater than the UV 
variability at A 1350. The EUV light curve shows contin- 
ual variations of a factor of ~ 2 on timescalcs of one or 
two days. The later EUVE light curve of Chiang et al. 
(2000) shows similar variations. For the most prominent 
event seen both in the EUV and UV spectra, Marshall 
et al. find a lag consistent with zero delay between the 
two bands. 

In Fig. 3 we show the correlation between the EUVE 
fluxes reported by Marshall et al. and the A1350 UV flux 
measured by the IUE and HST. We have averaged EUVE 
fluxes within ±ld of the UV flux measurement. From a 
least-squares fit we find that 

f E UV = 0-11/1350 - 97. (5) 

The range of EUV variability during the Marshall ct 
al. (1997) monitoring was a factor of two, while the UV 
range was only a factor of 40%. Equation (5) and Fig. 3 
imply that there is another component of UV flux that 
will remain even when the EUV flux goes to zero, at least 
on the timescale of the Marshall et al. observations. 



Support for a correlation of the EUV with longer wave- 
lengths is also provided by Dietrich et al. (2001), who 
made optical observations simultaneous with the Chiang 
et al. (2000) EUV and hard X-ray observations. Over the 
30-day monitoring period, the optical and EUV fluxes are 
correlated (see Fig. 4 of Dietrich et al.), but the EUV 
flux varied by a factor of > 6 while the observed opti- 
cal flux variations were only 20%. Dietrich et al. find 
optical-EUV lags consistent with zero, but with fairly 
large errors of ± 2 - 3 d . 

We used the UV flux rather than the optical flux to es- 
timate the what the observed EUV/soft X-ray flux would 
have been during the period we are considering because 
the difference in amplitude between the UV and EUV is 
not as great as between the optical and the EUV, and 
because the UV is closer to the EUV in log(v) than the 
optical is, and thus the UV has a higher probability of 
being physically connected to the EUV. The estimated 
observed EUV flux is not the true EUV flux because of 
the atomic absorption (see Steenbrugge et al. 2005 and 
references therein). We can get a better estimate of the 
true EUV and soft X-ray flux by scaling to the higher en- 
ergy side of this region, since the absorption problems are 
less. We therefore used the Walter et al. (1994) soft X-ray 
spectral shape relative to the comparatively absorption- 
free harder X-ray flux, which has been estimated above. 

2.5.3. Relationship Between the Soft and Hard X-ray 
Components 

The relationship between the soft and hard X-ray com- 
ponents provides an additional consistency check on our 
estimate of the soft component. As noted above, the hard 
X-ray variations lag the variations of the soft excess by a 
fraction of a day. From simultaneous ASCA and EUVE 
monitoring in 1996, Haba et al. (2003) found that, from 
day to day, the EUV and 2-10 keV X-ray fluxes were well 
correlated. The EUVE amplitude over a 10-day period 
was only 40% greater than the hard X-ray amplitude. As 
noted already, the 1998 simultaneous ASCA and EUVE 
monitoring by Chiang et al. (2000) showed a good corre- 
lation between the EUV and 2-10 keV X-ray fluxes, with 
the EUV amplitude being only 15 - 20% greater than the 
X-rays. Haba et al. did find, however, that on a shorter 
timescale (6 hours), there was a strong 40% dip in the 
EUV that had no counterpart in the hard X-rays (see 
their Fig. 2), but such rapid variability is not relevant 
for the emission lines we study here. 

Walter & Courvoisier (1990) fit the lowest energy ob- 
served with the thin Lexan filter on EXOSAT with a 
power law extrapolated from the hard X-rays and a soft- 
component. The soft-excess generally makes up 90% of 
the flux in the filter. In Fig. 4 we show the relationship 
between the soft-excess and the hard X-ray component. 
There is a similar relationship between the hard X-rays 
and the soft excess as between the UV and the soft ex- 
cess (see Fig. 3). Given that the 2-10 keV flux is corre- 
lated with the optical and UV on the BLR light-crossing 
timescale (see section 2.4), the amount of EUV variabil- 
ity implied by the relationship between the hard and soft 
X-ray fluxes is similar to that implied by the optical and 
UV.' 

3. THE GALACTIC REDDENING CORRECTION 
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Fig. 4. — Soft Excess (in counts) vs. Hard X-ray component (in 
units of 10 — 11 ergs s~ 1 cm -2 ) 

The observed continuum fluxes must first be corrected 
for Galactic extinction. Schlegel et al. (1998) used mea- 
surements of dust emission at 100 and 240 /im to cre- 
ate Galactic extinction maps, which have an accuracy of 
16%. From these the extinction, E(B-V), is 0.020±0.003 
mag for the direction of NGC 5548. From 21-cm observa- 
tions Murphy et al. (1996) have determined the neutral 
hydrogen column density, tih , towards NGC 5548 to be 
1.62 x 10 20 . If we adopt a standard dust-to-gas ratio of 
E(B - V) = 1.72 x 10" 22 n ff (Bohlin, Savage & Drake 
1978), we get E(B-V) = 0.028. If adopt a minimal dust- 
to-gas ratio of E(B - V) = 1.2 x 10" 22 n H (see Fig. 1 of 
Congiu et al. 2005) we get E(B-V) = 0.020. We can thus 
be confident that the Galactic reddening is relatively low 
and we adopt E(B -V) = 0.024 ± 0.004 

The Galactic extinction curve is well-determined for 
the near-IR, optical, and UV. For the far and extreme 
UV and soft X-ray regions the shape of the extinc- 
tion curve is not observationally determined so we used 
the interstellar dust model of Drainc (2003b) for R = 
Ay/E(B — V) = 3.1, which is a good fit to longer wave- 
lengths (Draine 2003a). 

4. HOST GALAXY CORRECTION 

The observed spectral energy distribution we have 
compiled is not the spectral energy distribution of the 
AGN alone: emission from the host galaxy is also in- 
cluded. Schmitt et al. (1997) have made a compilation of 
spectral energy distributions for normal galaxies, which 
enables us to subtract off the host galaxy contributions at 
various wavelengths. Romanishin et al. (1995) find that 
the galaxy component of NGC 5548 is a normal early- 
type spiral, with no conclusive evidence of star forma- 
tion in the bulge. We have therefore adopted the spiral 
galaxy SED of Schmitt et al. (1997). We scaled the spi- 
ral galaxy SED to a flux at 5200 A of 3.6 x 10" 15 ergs 
s" 1 cm" 2 A" 1 = 3.3 x 10" 26 ergs s" 1 cm" 2 Hz" 1 in the 
standard spectroscopic aperture of the I AW. 

5. INTERNAL REDDENING CORRECTION 

Gaskell & Benker (2007) have obtained extinction 
curves for individual AGNs by assuming that the unred- 
dened continua are the same. For most of the AGNs they 
consider this gives extinction curves which resemble the 
reddening curve in the local solar neighborhood minus 
the A2175 amorphous carbon feature, but which are flat- 
ter in the UV (see also Czerny et al. 2004 and Gaskell 
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Fig. 5. — The extinction curve for NGC 5548 (filled circles) 
compared with the mean AGN reddening curve of Gaskell & Benker 
(2007) (solid curve) and an SMC reddening curve (dashed line). 

et al. 2004). For one of the AGNs they consider they 
find a steeply-rising, SMC-type reddening curve. We ob- 
tained the extinction curve for NGC 5548 in the same 
manner as Gaskell & Benker by comparing our observed 
spectrum of NGC 5548 with the observed blue spectrum 
of the relatively unreddened AGN 3C 273. As can be 
seen in Fig. 5, the resulting reddening curve is consistent 
with the mean AGN reddening curve of Gaskell & Benker 
(2007) but is inconsistent with an SMC reddening curve. 
Because of the consistency with the standard AGN red- 
dening curve we have adopted the latter for NGC 5548. 
The normalization of NGC 5548 to the standard AGN 
curve gives E(B — V) = 0.17 mag. This has a formal 
uncertainty of ~ ±0.01 mag., but the real uncertainty is 
probably several times this because of the uncertainty in 
the host galaxy correction for 3C 273. 

If we assume the standard Galactic dust-to-gas ratio 
of Bohlin et al. (1978) within NGC 5548, E(B-V) = 0.17 
corresponds to a total hydrogen column of 10 21 cm" 2 . 
This is the same as the median column density Walter & 
Courvoisier (1990) obtain. Steenbrugge et al. (2005) get 
total best-fit column densities for their warm absorber 
components of 3 x 10 21 cm" 2 which would imply E(B- 
V) ~ 0.5 with a standard gas/dust ratio but this is not 
necessarily inconsistent with our lower reddening since 
not all warm-absorber gas will have dust associated with 
it. 

The Gaskell & Benker (2007) reddening curve only ex- 
tends to Lyman a. They find that the slope of AGN spec- 
tra in the far UV (A < Ly a) is not well correlated with 
E(B-V) determined from UV-optical spectra (see their 
Fig. 7) and they propose that the scatter in the far-UV 
slope is the result of a small amount of additional redden- 
ing (E(B-V) ~ 0.03 mag.) by a SMC-like dust since this 
spectral region is very sensitive to such dust. The bluest 
AGNs have a spectral index from ~ 1200 A to ~ 800 A of 
~ 0.75. We have adopted this as the unreddened spectra 
slope in this spectral region. 

The mean AGN extinction curve is uncertain at high 
energies. We therefore extrapolated the extinction curve 
to the X-ray region using the average of the Milky Way 
R = 3.1 and R = 5.5 curves calculated by B. Draine. 
These is based on the carbonaceous - silicate dust model 
of Weingartner & Draine (2001) and Li & Draine (2001). 
The graphite and silicate optical properties are given in 
Draine (2003a, b), and include structure around X-ray 
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Fig. 6. — Estimated spectral energy distribution of NGC 5548 
around JD 2449120. The upper dashed line represents the raw ob- 
served or estimated continuum before any reddening or host galaxy 
corrections. The thick solid line and the dots show the corrected 
observed continuum from the optical to the X-ray region including 
the interpolations in the EUV gap. This represents the SED as 
it would be seen from the earth with no extinction and no host 
galaxy light. The straight line shows the gap. The IR observations 
are shown as dots, while the lower dashed lines show the IR at 7 
Id (thick dashed) and at 80 Id (thin dashed). The vertical scale is 
arbitrary. 



TABLE 2 
De-reddened Continuum Fluxes 





Log(*AF„) a 


Log(^) 


Log(i^) a 


12.5 


-10.19 


16.5 


-9.79 


13.0 


-10.07 


17.0 


-10.47 


13.5 


-10.17 


17.5 


-10.63 


14.0 


-10.20 


18.0 


-10.81 


14.5 


-10.25 


18.5 


-10.51 


15.0 


-9.87 


19.0 


-10.52 


15.5 


-9.51 


19.5 


-10.74 


16.0 


-9.38 


20.0 


-11.79 



a vFv in ergs s 1 cm 2 at the earth 



absorption edges. 

6. ADOPTED SPECTRUM 

The only remaining uncertainty in our derived spec- 
trum (see Fig. 6) is the "EUV gap" . After our reddening 
corrections the v¥v spectrum is rising on either side of 
the gap. We interpolated across the gap by fitting a 6th- 
order polynomial to the points on either side of the gap. 
This gives a good fit on either side of the gap and we 
consider this to be the most likely interpolation through 
the gap. We also did a simple linear interpolation across 
the gap and this is also shown in Fig. 6. The latter is 
most unlikely to be the correct interpolation, but this ex- 
treme case shows the maximum uncertainties introduced 
by our extrapolation. For convenience we also tabulate 
our final continuum in Table 2. This gives the continuum 
as it would be seen from the earth with no extinction and 
no host galaxy light contamination. 



7. EMISSION-LINE FLUXES 

In Table 3 we give the observed emission-line fluxes for 
the highest state during the period we consider. The 



TABLE 3 

Emission-Line Fluxes 



Line 




Obsd. 


MW de-red 


Int. de-red 


Lya 


A1216 


480 


600 


2300 


N V 


A 1240 


93 


120 


450 


O I 


A1304 


30 


40 


140 


O IV] + Si IV 


A1400 


90 


110 


430 


C IV 


A1549 


730 


870 


3200 


He II 


A1640 


75 


89 


320 


C III] 


A1909 


135 


160 


550 


Mg II 


A2798 


130 


150 


400 


He II 


A4686 


15 


16 


30 


H/3 


A4861 


62 


70 


120 


He I 


A5876 


11 


12 


18 


Ha 


A6563 


250 


270 


390 



Ly a, C IV, and CIII] fluxes are from Korista et al. 
(1995), the H/3 flux is from Wanders & Peterson (1996) 6 , 
and the He II A1640 flux is from Bottorff et al. (2002). 
As Bottorff et al. discuss, there is a large uncertainty 
(~ 30%) in this flux. There are also even larger un- 
certainties in measuring He II A4686 (see Appendix in 
Bottorff et al. 2002) and, in addition, the He II A4686 
line was not measured in the 1993 campaign. The value 
given in Table 1 was scaled from the 1989 campaign re- 
sults. Ha was also not observed in the 1993 campaign 
so the strength given in Table 1 comes from assuming 
that Ha/H/3 w 4.1 (Shapavalova et al. 2004). The slight 
variations in this ratio with time (see Shapavalov et al.) 
are unimportant for our analysis. He I A 5876 was also 
not observed during the 1993 campaign, so in Table 3 
we have simply scaled it from H/3 using the median He 
I A5876/H/3 ratio from Osterbrock (1977). As can be 
seen from Shapavalova et al. (2004) He I A5876 varies 
substantially. Table 3 also gives the fluxes corrected for 
Galactic and internal reddening as discussed above. 

8. REDDENING FROM LINES 

Emission-line ratios can potentially give an indepen- 
dent check on reddening derived from the continuum. 
The He II recombination lines give the most intrinsically 
reliable line ratios, since they are mostly independent 
of physical conditions, but, as has been extensively dis- 
cussed by Bottorff et al. (2002), measuring the He II 
A1640/A6886 ratio is very difficult (see their Appendix), 
so while the He II fluxes in Table 3 are consistent with the 
adopted reddening, the He II line ratio does not provide 
a reliable constraint. 

The Balmer decrement is sensitive to physical condi- 
tions and varies both with projected velocity (see Sned- 
den & Gaskell 2004, 2007a,b)' and with time, but the 
relative integrated Ha/H/3 ratio is an approximate indi- 
cator of reddening (e.g., De Zotti & Gaskell 1985). Sha- 
pavalova et al. (2004) find Ha/H/3 varying from 3.5 to 
4.3, which is not too different from Osterbrock (1977). 
If we adopt an unreddened decrement of 3.0, then the 
NGC 5548 decrements give E(B-V) - 0.17 - 0.36, both 
which are consistent with the reddening we get from the 
continuum. 

9. PHOTOIONIZATION MODELLING 

6 See Footnote 1 
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After determination of the continuum we ran models 
using the photoionization code CLOUDY, version 96.01 
(Ferland et al. 1998) to determine the emission-line spec- 
trum. In order to match the source-to-cloud distances 
found from reverberation mapping (Korista et al. 1995) 
and our estimated 80 Id inner edge of the dusty torus, we 
ran CLOUDY models corresponding to source distances 
of 7, 14, 21, 30, 50, and 80 Id. The continua at each 
of these radii were calculated as described above. Since 
there will certainly be a range of densities in the BLR, 
we ran models with hydrogen densities, n#, of 10 95 , 
lO 10 5 , and 10 11 " 5 cm" 3 . Because they are located at 
the very centers of galaxies, the metalicities of AGNs 
are almost always supersolar (Hamann & Ferland 1999), 
so we adopted a metal abundance of three times solar. 
All models run were completely optically thick models 
and the calculations were continued far deeper into the 
clouds than any of the lines we are interested in were 
being produced. 

In addition to these models we ran additional models 
to investigate the effect of not reducing the IR emission, 
and of varying the metal abundances from solar to six 
times solar. We also investigated the effect of significant 
(1000 km s" 1 ) (non-dissipative) internal turbulence (see 
Bottorff et al. 2000 and Bottorff & Ferland 2002). 

10. BROAD-LINE COVERING FACTORS 

In Table 4 we present the covering factors needed to 
produce the observed line fluxes. In order to show the 
effect of the radial distribution of the gas for each line we 
give the covering factors needed as a function of radius 
under the assumption that each line is only produced in 
gas at this radius. We show the effects of varying den- 
sity, source distance, metal abundance, and continuum 
shape, and of adding turbulence. As well as giving cov- 
ering factors using our standard continuum ( "adopted" ) 
we also show the effects of using the simple linear inter- 
polation across the EUV gap ("linear"), of using just the 
observed SED ("observed") with no corrections at all, 
and of not allowing for the BLR being far from the IR 
emission ("Extra IR"). As is to be expected, changing 
the EUV continuum shape has the biggest effect on the 
He II recombination lines (see, for example, Mathews & 
Ferland 1987). Adding extra IR heating increases the 
strengths of all collisionally-excited lines. 

Looking at the density and radius dependencies of the 
covering factors for our standard model it can be seen 
that, as is well known (see, for example, Laor 2007), 
the production efficiency of most lines decreases at high 
ionizing flux levels. The efficiency generally increases 
with density for many lines except where lowering the 
resultant lowering of the ionization reduces the line flux 
(e.g., N V), or a line is collisionally suppressed (e.g., C 
III]). 

The observed integrated line fluxes of an AGN are the 
sums of the contributions from all clouds at all densities 
and all radii (cf. the "LOC" model of Baldwin et al. 1995 
and Korista & Goad 2000). In the bottom line of Table 
4 we give the covering factors for an integrated model 
based on our standard continuum model, with a uniform 
range of densities from 10 95 to 10 115 cm" 3 , and with 
uniform covering factor as a function of radius. There is 
generally good agreement in the resulting covering fac- 
tors which range from 17-56% (ignoring He I A5876 for 



which the observed flux is highly uncertain) and with 
an average is 40%. This average covering factor can be 
lowered somewhat by weighting the density distribution 
towards higher densities, and by having the covering fac- 
tor increase with radius, but it can be seen from Table 
4 that an average integrated covering factor of < 20% 
is not attainable. The only way of lowering the aver- 
age covering factor to < 20% is by including substantial 
internal turbulence. Further lowering of the covering fac- 
tor can be achieved by making the turbulence dissipative 
to provide additional heating for the low-ionization lines 
(Bottorff & Ferland 2002). 

11. THE TORUS COVERING FACTOR 

The dust covering factor, £ldust, can be found from the 
ratio of the bolometric luminosity, Li> i, to the luminos- 
ity, in the mid- infrared, Lj/j. Since we have the complete 
spectrum for NGC 5548 we can perform a direct integra- 
tion to obtain Lb i- We show the results for different 
continuum assumptions in Table 5. The observed case 
ignores reddening and the host galaxy correction. The 
"low big bump" case is the simple linear interpolation 
across the EUV gap and the "high big bump" case is 
the more likely interpolation (see Fig. 6). The covering 
factors given in Table 5 assume that the torus has an 
unobstructed view of the higher-energy continuum. We 
will argue below that this is probably not the case. 



12. VARIABILITY TRANSFER FUNCTIONS 

The radial distribution of material in the BLR can only 
be inferred indirectly from reverberation mapping (see 
Peterson 2006 for a recent review). The most widely used 
measure of the size of the BLRs of AGNs is the "lag" 
given by the centroid of the cross-correlation function 
(Gaskell & Sparke 1986, Gaskell & Peterson 1987), but 
this has always been known to be biased towards gas at 
small radii. A better indication of the distribution of the 
gas is given by the transfer function \P(t), which is the 
response of the system at a delay t to a delta function at 
r = 0. A detailed comparison of observed and theoretical 
transfer functions is beyond the scope of this paper and 
will be given elsewhere, but in Figs. 7 - 9 we compare 
observed transfer functions of C IV and H/3 with simple 
theoretical transfer functions based on the parameters of 
our photoionization modelling. 

Recovering the ^(r) from observed continuum and 
line-flux time scries is difficult. In addition to uncer- 
tainties arising from noise in the data (which are most 
significant for weak lines like those of He II), the derived 
^(r) depends on the sampling, the method used, and as- 
sumptions that always need to be made to recover \I/(t) 
in practice. Assumptions need to be made about the 
smoothness of ^(t) and boundary conditions including 
the contribution of relatively non-varying components. 
It also needs to be recognized that some methods (e.g., 
the maximum-entropy method) only allow ^(r) to be 
positive, so noise in ^(t) gives spurious positive con- 
tributions at large delays. Since AGN variations arc 
commonly quasi-cyclical on the timescale of almost all 
monitoring campaigns, aliasing at multiples of the quasi- 
period of the continuum variations will also make spu- 
rious contributions to ^(r) at large delays. Because 
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TABLE 4 

Covering Factors 



Model 


log(n w ) 


r(ld) 


Lya 


N V 


O I 


A1400 


C IV 


He II 


C III] 


Mg II 


H/3 


He I 


Adopted SED 


9.5 


7 


74% 


157% 




192% 


223% 


25% 


77% 


212% 


426% 






9.5 


14 


28% 


86% 


417% 


113% 


131% 


17% 


32% 


306% 


80% 


123% 




9.5 


21 


17% 


60% 


196% 


85% 


106% 


19% 


21% 


119% 


44% 


66% 




9.5 


30 


13% 


43% 


135% 


70% 


83% 


21% 


16% 


68% 


28% 


43% 




9.5 


50 


10% 


28% 


106% 


54% 


55% 


22% 


11% 


28% 


18% 


29% 




9.5 


80 


9% 


21% 


87% 


44% 


35% 


24% 


10% 


15% 


15% 


30% 




10.5 


7 


35% 


43% 


173% 


77% 


145% 


15% 


81% 


526% 


218% 
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Observed 


10.5 


21 


64% 


55% 


147% 


69% 


71% 


100% 


77% 


51% 


115% 


107% 


Adopted 


10.5 


21 


15% 


22% 


31% 


64% 


70% 


22% 


50% 


101% 


57% 


50% 


Linear 


10.5 


21 
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31% 
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84% 


74% 


41% 


75% 


100% 


65% 


62% 


Adopted 


10.5 


7 


35% 


43% 


173% 


77% 


145% 


15% 


81% 


526% 


218% 


181% 


Extra IR 


10.5 
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908% 
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3 solar 
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6 solar 
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3 solar 
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14 
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101% 
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200% 
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6 solar 
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10.5 


21 


15% 
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31% 


64% 


70% 
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50% 


101% 
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50% 


Turb. 


10.5 


21 


8% 


16% 


58% 


20% 


27% 


21% 


26% 


12% 


11% 


25% 


LOC Model 






17% 


34% 


61% 


51% 


56% 


22% 


37% 


28% 


36% 


13% 



TABLE 5 

ESTIMATED DUST COVERING FACTORS 



Continuum Shape &dust 

Observed (uncorrected) 61% 

Internally de-reddened with low big bump 25% 

Internally de-reddened with high big bump 19% 



of these mcthodolical differences, a variety of different 
transfer functions have been published from the same 
data sets even for strong lines. 

As is well known, ^(r) for a thin spherical shell of ra- 
dius r is a boxcar from r = to r = 2r/c. The transfer 
functions for other spherically-symmetric distributions 
can be constructed by adding such boxcar functions. For 
example, a filled sphere has ^(t) declining monotonically 
from a peak at r — 0. A hollow sphere of inner radius rj„ 
has ^(t) constant until r = 2r in /c. If the clouds emit 
isotropically, the only way to avoid 'J'(t) having a maxi- 
mum at r = is to have a non-spherical geometry. Most 
type-1 AGNs are observed close to pole on so we mod- 
elled non-sphcrically symmetric distributions as spheres 



with polar cones cut out of them viewed from on axis. A 
single thin sphere with a cap cut out around each of the 
poles gives \£(r) as a box car of width ± sin(#)r/c cen- 
tered on r/c, where 9 is the maximum angle the clouds 
subtended above the equatorial plane, ^(t) for radially- 
thick distributions is simply the sum of these box cars. 
We calculated simple theoretical AGN line transfer func- 
tions, smoothed to the observed resolution, for clouds in 
a distribution which is truncated by cones on the symme- 
try axis, but is otherwise spherically symmetric, ^(r) for 
geometries where the observer is slightly off axis, as will 
be the case for most AGNs, are mostly indistinguishable 
from on-axis cases. Off-axis cases can be divided into two 
classes: when the observer is inside the cone, and when 
he or she is outside the cone. If the opening of the cone is 
wide, the former case will be indistinguishable from the 
on-axis case, and if the opening of the cone is narrow, 
the latter case will be indistinguishable from spherical 
symmetry. 

Our models have only three main parameters: the half- 
opening angle of the polar cone, 90deg— 9, the outer ra- 
dius, and the radial variation of responsivity which we 
parameterized as a simple power laws. The radial vari- 
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Fig. 7. — Theoretical (solid) and observed (dashed) transfer 
functions for H/3 in NGC 5548. The observed transfer function is 
from Home, Welsh, & Peterson (1991). 

ation of the responsivity depends on three things: the 
response of a particular line to changes in the ionizing 
continuum, the emissivity per cloud as a function of ra- 
dius, and the dependence of the covering factor (number 
of clouds) as a function of radius. The functional form 
of the response of the flux in a given line from a cloud to 
changes in the continuum is expected to remain nearly 
constant over the range of conditions considered. The ra- 
dial variation of responsivity thus depends on the change 
in emissivity, which is given by our photoionization mod- 
els, and the unknown dependence of the covering factor 
on radius. For simplicity we assumed that the covering 
factor, which has to be the same for all lines coming from 
a given type of cloud, has no radial gradient, and that 
the radial variation in responsivity thus depends solely 
on the radial variation of emissivity given by the pho- 
toionization models. The remaining parameters for each 
transfer function are the opening angle and the outer ra- 
dius. Although these assumptions are very simple, they 
nonetheless give useful insights into the behavior of trans- 
fer functions. 

12.1. The HP Transfer Function 

The photoionization models show that the emissivity 
of H/3 rises as ~ r +1 . Home, Welsh, Peterson (1991) 
found that ^(t)hp (which we show as a dotted line in 
Fig. 7) does not peak at zero delay. Cacklett & Home 
(2006) show that over the course of the 13-year IAW 
monitoring of NGC 5548 $(t)h0 maintained the same 
shape regardless of the mean flux level of NGC 5548. As 
outlined above, such a transfer function demands a pre- 
dominantly non-spherical distribution 7 seen from within 
the polar opening. As can be seen in Fig. 7, we can 
readily reproduce the shape of ^{t)h/3 with our simple 
models. The symmetry of ^(t)h0 requires 9 < 25deg. 
This corresponds to a covering factor of 40% or less, 
and thus is in good agreement with the results in Ta- 
ble 4. For such a distribution, the position of the peak 
of ^(t) depends only on the outer radius of the distribu- 

7 Ferland et al. 1992 show how a spherical distribution can pro- 
duce the H/3 transfer function if the emission is preferentially back 
towards the source, but this is unlikely since the Ly a transfer 
function in NGC 5548 is incompatible with such asymmetric emis- 
sion. Also, the double-humped H/3 profile is independent evidence 
against spherical symmetry. 
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Fig. 8. — Theoretical (solid) and observed (dashed) transfer 
functions for He II A1640 in NGC 5548. The observed transfer 
function is from Krolik et al. (1991). 

tion. The outer edge of the BLR, r out is expected to be 
the dust sublimation radius (Netzer & Laor 1993, Laor 
2007) which is quite well defined. From consideration 
of IR reverberation mapping in section 2.2 we estimated 
that rdust ~ 80 Id, but we can get a more precise esti- 
mate of r ou t by fitting ^>(t)h/3- At the optical flux level 
of NGC 5548 during the high state we are considering, 
Cackett & Home (2006) find the peak of the ^(t) H to 
be at 19 Id, with lower and upper quartile uncertainties 
of -2 and +5 Id respectively. This is consistent with the 
lag determined by the Korista et al. (1995). Fitting our 
model transfer function gives r ou t = 60 Id with the same 
fractional uncertainty as in the position of the observed 
peak of ^(t)h/3- 

12.2. The He II Transfer Function 

The transfer function of He II was determined by Kro- 
lik et al. (1991). It is not as well determined as that of 
H/3, but, as can be seen in Fig. 8, it shows a strong peak 
at zero delay. The emissivity of He II, unlike that of H/3, 
declines slowly with radius as <~ r~ 0A . In our models this 
declining emissivity automatically gives ^(r) peaking to- 
wards zero delay. The radical difference in the shape of 
^(t)h/3 and ^(r)Heii thus primarily arises from pho- 
toionization physics (how clouds emit as a function of 
incident flux) rather than a profoundly different spatial 
distribution. Nonetheless, there must be some difference 
in the geometries of the clouds dominating the H/3 emis- 
sion and those dominating the He II emission. If we 
are to keep the same flattened distribution derived from 
the ^>(t)h/3 then ^{r)Heii declines too slowly with in- 
creasing delay. To make our model match the observed 
^(T)Heii requires us to reduce r out to ~ 30 days. How- 
ever, if the clouds are directly exposed to the ionizing 
continuum (as in our CLOUDY models), we still have 
He II emission wherever there is H/3 emission, and the 
He II emitting cloud distribution would not be truncated 
more than the H/3 emitting cloud distribution. In this 
case another way of matching the shape of the He II 
transfer function (see Fig. 8) is to increase the covering 
factor. For pole-on viewing a covering factor of over 90% 
is needed (9 > 60deg). A more likely possibility is that 
the covering factor is only slightly greater than for H/3 
(and hence in agreement with the results of Table 4), but 
we are now outside the hollow polar cone. 
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Fig. 9. — Theoretical and observed transfer functions for H/3 in 
NGC 5548. The theoretical transfer function (see text) is shown 
as a solid line. The squares are the transfer function derived with 
the SOLA method from the 1993 campaign (Wanders et al. 1995) 
and the dashed line is the transfer function derived from the 1989- 
90 campaign by Home et al. (1991) using the maximum entropy 
method. 

12.3. The C IV Covering Factor 

In Fig. 9 we show &(t)civ estimated at two epochs 
by two different methods by Wanders et al. (1995) and 
Home et al. (1991). Additional estimated C IV transfer 
functions can be found in these papers and also in Done & 
Krolik (1996). There is general agreement in ^{t)civ f° r 
r < 20 d for both epochs, but more uncertainty for r > 
30 d. For C IV the radial dependence of the emissivity is 
intermediate between that of H/3 and He II. In the model 
shown in Fig. 9 we adopt a flat radial dependence. If 
we again assume r out <~ 60d, as estimated from ^(r)Hfi, 
then to fit the observed ^(t)civ we need a high covering 
factor as for He II. Again this could really be a lower 
covering factor with us looking through the equatorial 
clouds. This is quite likely to be the case for NGC 5548 
since the appearance of a double-humped, disk-like H/3 
profile implies that we are not viewing it pole-on. It is 
also well known (see Arav et al. 2003, Steenbrugge et 
al. 2003 Steenbrugge et al. 2005, and refs. therein) that 
there is C IV absorption in NGC 5548. 

13. DISCUSSION 
13.1. Is There an Energy-Budget Problem? 

We have considered in detail the overall SED of 
NGC 5548 during the high state at JD 2449120 and 
argue that although not all spectral regions were ob- 
served, we can be fairly confident of the overall SED. The 
main assumption still to be confirmed is the suggestion 
of Gaskell & Benker (2007) that the apparent turndown 
in the SED below Ly a is not intrinsic but results from 
a small amount of SMC-like dust. If this is incorrect, it 
only has a slight effect on our estimated covering factors. 

The He II recombination lines, which are a very good 
indicator of covering factor since they are very insensi- 
tive to physical conditions, only require a covering fac- 
tor of ~ 25% over a wide range of conditions with our 
adopted model. Other lines depend on physical condi- 
tions more but in our integrated ("LOC") model only 
two lines, O I A1304 and C IV require covering factors 
of more than 50%. Only the He II flux depends signifi- 
cantly on our assumption that there is a small amount of 



SMC-like reddening. O I A 1304 is very sensitive to the 
mctalicity and a slight increase in metalicity is sufficient 
to bring the O I covering factor in line with the other 
lines. High metalicity also lowers the covering factors 
needed for N V, the O IV] + Si IV A1400 blend, and, to 
a lesser extent, C III]. Increasing the metalicity does not 
help reduce the C IV covering factor since it is a major 
coolant. What does help reduce the C IV covering fac- 
tor most is an increased contribution of higher density 
gas. Higher density gas generally increases the emission 
efficiency of most lines. 

From Table 4 it can be seen that with what we con- 
sider to be the best continuum, most of the lines can 
be explained by a covering factor of ~ 40% or some- 
what less. This is not significantly lower than covering 
factors obtained in the past with weaker EUV continua 
(see, for example, Shields & Ferland 1993; Korista et 
al. 1998; Goad & Koratkar 1998; Kaspi & Netzer 1999; 
Korista & Goad 2000), so our main point is that even 
though we have argued that because of reddening there 
are more photons in the EUV continuum than hitherto 
recognized, is it not possible to lower the covering factor 
significantly. In particular, it is not possible to lower it 
to the level required by the Lyman limit absorption con- 
straints. The reason why reddening does not solve the 
problem is that the lines have to be de-reddened as well 
as the continuum. 

Independent support for the BLR covering factor we 
derive comes from the Fe Ka line equivalent width. Chi- 
ang et al. (2000) derive a Fe Ka covering factor of 20-25% 
for NGC 5548. 

While all lines can be explained with covering factors 
significantly less than 100%, and so there is formally no 
"energy-budget" problem, we are unable to lower the 
necessary covering factors by an order of magnitude. In 
the old scenario of clouds uniformly covering the cen- 
tral source, a covering factor of ~ 40% (or even 20%) 
would still be a problem because Lyman continuum ab- 
sorption is certainly not detected in 40% of AGNs. We 
believe that the solution to this problem, as proposed 
by Maiolino et al. (2001c), is that the clouds do not 
cover the source uniformly. Instead, they are in an ax- 
imuthally symmetric distribution and we are looking in 
through a hole at the pole. Important support for this 
comes from the line- continuum variability transfer func- 
tions discussed above. We have shown that the geometry 
implied by the well-measured H/3 transfer function is in 
good excellent quantitative agreement with this scenario. 
For other lines (e.g., C IV and He II) a less flattened dis- 
tribution is needed. Additional support for a flattened 
geometry comes from the alignment of the polarization 
vector parallel to the axis of symmetry (see discussion in 
Goosmann & Gaskell 2007), and from polarization rever- 
beration (Gaskell et al. 2007). 

Since our deduced covering factors can be reconciled 
with line transfer functions for near pole-on viewing an- 
gles, we do not believe that there is an energy-budget 
problem for NGC 5548. This also makes some previously 
suggested solutions to the problem unnecessary. Given a 
flattened geometry there is no need for additional heating 
(Dumont et al. 1998) or for a large amount of turbulence 
(Bottorff et al. 2000, Bottorff & Ferland 2002) to reduce 
the covering factor, although such factors could still be 
relevant. 
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Fig. 10. — Schematic representation of BLR geometries, (a) Tra- 
ditional model where BLR clouds are distributed quasi-sphcrically 
and mostly do not shield distant clouds or the torus, (b) Shielding 
model where the clouds are concentrated in the equatorial plane. 
The torus and distant clouds see the central source shielded by 
clouds closer to the center. (In both sketches the BLR clouds are 
not shown to scale and the number of clouds has been reduced by 
many orders of magnitude.) 

13.2. A New BLR Model 

If the majority of the BLR gas is located within 
±25 deg of the equatorial plane, and the overall cover- 
ing factor f2/47r w 0.4, then photons leaving the central 
source in a direction close to the equatorial plane have 
a close to 100% probability of striking BLR gas. By 
the time one gets to the inner edge of the dusty torus 
the spectrum has been substantially modified by passing 
through the inner BLR. Such modification of the spec- 
trum has already been considered by Holt et al. (1980) 
and modelled for NGC 4151 by Ferland & Mushotzky 
(1982). The effect of shielding is similar to taking a 
standard cloud, cutting it in half in a plane perpendicu- 
lar to the direction of the incoming photons and moving 
the outer part of the cloud further away. The BLR sce- 
nario we require resembles the geometry considered by 
MacAlpine (1972). In Fig. 10 we contrast our proposed 
shielded BLR geometry with the widely-used geometry 
where shielding is assumed to be negligible. 

As discussed in 12.2 and 12.3, the differences in the 
transfer functions of C IV and H/3 require either that 
the high-ionization emission comes from a smaller region, 
or that it has a higher covering factor, or both. It has 
in fact long been argued that there are effectively two 
separate types of cloud: high-ionization clouds and low- 
ionization clouds (Gaskell 1987; Collin-Souffrin & Lasota 
1988; Gaskell 2000). Not only are differences in radii 
implied by reverberation mapping, but the higher ve- 
locity dispersions of high-ionization lines (Osterbrock & 
Shuder 1982; Shuder 1982) show that the high-ionization 
gas is substantially closer (see Peterson & Wandel 1999). 
The blueshifting of high-ionization lines (Gaskell 1982) 



TABLE 6 

MEASURED AND PREDICTED LAGS FOR NGC 
5548 



Line 




Observed^ 


LOC b 


Shielded 


Ly a 


A1216 


12 


8.5 


12 


N V 


A1240 


4 C 


3.8 


2.4 


O I 


A1304 


40 c 




43 


Si IV + O IV] 


A1400 


12 


5.8 


15 


C IV 


A1549 


8 


8.0 


8 


He II 


A1640 


3 d 


7.3 d 


3.3 


C III] 


A1909 


35 


11.5 


28 


Mg II 


A2798 


> 34 


16 


28 


He II 


A4686 


3 C 


7.3 d 


3.4 


H/3 


A4861 


20 f 




18 


He I 


A5876 


13 g 




16 


Ha 


A6563 


22.5 g 




19 



a From Clavel et al. (1991) except as noted. 

b Based on Korista & Goad (2000) except as noted. 

c From Krolik et al. (1991). 

d From Bottorff et al. (2002). 

e Set equal to A1640 lag (see footnote 8). 

f From Peterson et al. (1991). 

s Scaled from H/3 using square of mean ratio of FWHMs 
from Osterbrock & Shuder (1982). 

requires that the gas producing the high-ionization lines 
is spatially separated from the low- ionization gas. Fur- 
thermore, for NGC 5548 and other objects we find that 
the velocity dependence of physical conditions is different 
for the high- and low-ionization gas (Snedden & Gaskell 
2004; Nazarova et al al. 2004, 2006; Snedden & Gaskell 
2007a,b). 

There have been attempts (notably by Korista & Goad 
2000) to model this stratification in detail with a con- 
tinuous distribution of cloud properties as in the LOC 
model, but there is a major difficulty in quantitatively 
matching the range of observed lags and line widths. As 
Korista & Goad (2000) point out, there are problems at 
both the high- and low-ionization ends: He II and N V 
have much smaller lags than predicted, while C III] and 
Mg II have significantly longer lags than predicted. He II 
presents another problem. The LOC models (Bottorff et 
al. 2002) predict that the radius of the He II emitting re- 
gion should be only 10% smaller than the C IV emitting 
region, but as Bottorff et al. point out, the observations 
for NGC 5548 (in both 1989 and 1993), and for NGC 3783 
(Reichert et al. 1994), Fairall 9 (Rodriguez-Pascual et al. 
1997), NGC 7469 (Wanders et al. 1997), and 3C 390.3 
(O'Brien et al. 1998) all give a lag of He II A1640 three 
or four times smaller than the C IV lag. 8 In Table 6 we 
compare the observed lags for NGC 5548 (normalized to 
a C IV lag of 8 days) with the predictions of the LOC 
model predictions. The LOC predicted lags have also 
been scaled to a C IV lag of 8 days. It can be seen the 
LOC model only predicts a factor of four variation in lag 
while the observations give three times that. 



The shielding model shown in Fig. 10b naturally ex- 

8 Bottorff et al. (2002) also point out a problem with He II 
A4686 appearing to give a lag longer than that of He II A1640. 
The explanation of the longer He II A4686 delay is probably that 
a small amount of contamination by Fe II and the blue wing of H/3 
makes the He II A4686 lag appear longer. See Gaskell (2007). 
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plains the increased stratification. The relative order of 
ionization stages in a photoionized cloud is very insensi- 
tive to changes in the density or the ionizing flux; only 
the scale factor changes. We can therefore approximate 
a shielded distribution of clouds by taking a single cloud 
in CLOUDY and pretending that it has been broken into 
smaller pieces spread uniformly throughout the BLR. We 
simply linearly transform the depths into the cloud from 
the side facing the ionizing source into radii away from 
the central source. For different lines of sight and for dif- 
ferent densities, cloud sizes, and filling factors the scale 
factor will be different, but the order of stages of ion- 
ization remains the same. The weighted mean radii of 
maximum emission of each line will therefore remain in 
the same ratios as the depths into a single cloud. One 
important difference from a single cloud is that photons 
can now escape from the distribution much more easily. 
We treated this in the same manner as dissipationless 
turbulence. 

To quantitatively test the shielded BLR model, we took 
the 21 Id, 10 10 - 5 cm" 3 CLOUDY model with a 1000 km 
s -1 turbulence velocity, and scaled the radius of peak 
C IV emission to match the observed 8-day lag. The 
only other free parameter is the outer edge of the BLR 
which we took to be the 80 Id inner edge of the torus. We 
show our predicted lags in Table 6. As can be seen this 
give a better quantitative agreement with the observed 
lags (±20%) than the LOC model. 

13.3. The Torus Covering Factor 

Netzer & Laor (1993) proposed that the outer edge of 
the BLR is set by the dust-sublimation radius. Our mod- 
eling of transfer functions requires a fixed outer radius of 
~ 60d, and we have shown that this is in agreement with 
the radius of the hot dust measured by IR reverberation 
mapping. 

We consider the most likely dust covering factor to be 
19% (see section 11). This is lower than almost all the 
covering factors we deduced for the various emission lines 
in Table 4 and would imply that the dust has a flatter 
distribution than the gas. 

The dust covering factor we find for NGC 5548 is simi- 
lar to those of other AGNs. Maiolino et al. (2007) discuss 
the dust covering factors for a wide range of AGNs. They 
estimate the dust covering factor by assuming a ratio of 
LboI = 7z/L„(A5200). However, the study of Shang et 
al. (2005) implies that Lb i = 15^i J/ (A5200) is more ap- 
propriate, so that the covering factor can be calculated 
from 

O w 0.18^(6. 7^/^(51001) (6) 

For NGC 5548 this gives VL/Ait rts 0.18 which is compa- 
rable to dust covering factors Maiolino et al. (2007) find 
for AGNs of comparable luminosity if we correct for our 
updated scaling of ^L„(A5200) to Lb i 

Although the deduction of Sldust /47r w 0.2 is straight 
forward and consistent with the covering factors we de- 
duced for low-ionization lines, there is a problem: our 
NGC 5548 torus covering factor and similarly estimated 
covering factors for other AGNs are smaller than the 
torus covering factors implied by the observed ratio of 
type-l/type-2 AGNs. The ratio of type-l/type-2 AGNs 
in hard X-ray samples (see Fig. 7a of Maiolino et al. 2007) 



implies that at the luminosity of NGC 5548, Sldust/^K = 
0.6. A further problem of having fltorus < &blr is that 
there should be many more objects where we see strong 
low-ionization absorption by the BLR (Maiolino et al. 
2001c). While there are low-ionization broad absorption 
line AGNs, this absorption is not arising from the main 
BLR because it is blueshifted, and hence arising from 
outflowing gas. Because the variability of the blue wings 
of broad emission lines does not lead that of the red wings 
(Gaskell 1988; Koratkar & Gaskell 1989, Crenshaw & 
Blackwell 1990) the main BLR is not outflowing but is 
gravitationally bound. Observations indicate in fact that 
the BLR is inflowing slightly (Gaskell, 1988, Koratkar & 
Gaskell 1989, Koratkar & Gaskell 1991a,b, Korista et 
al. 1995; Done & Krolik 1995; Gaskell & Snedden 1997, 
Welsh et al. 2007). The problem of the lack of BLR ab- 
sorption at the systemic velocity requires Cltorus < lr 
We believe that our proposed BLR geometry offers the 
solution to these two problems: the BLR is almost cer- 
tainly co-planar with the torus, so the latter does not 
have an unobscured view of the central source. There- 
fore, the energy received by the torus is substantially re- 
duced, and hence the torus has to subtend a larger solid 
angle to explain the observed IR luminosity. If we assume 
that all of the energy from the Lyman limit to 1 keV is 
absorbed, the dust covering factor increases by a factor 
of 2.7 to 53%, which is consistent with the mean torus 
covering factor expected from the ratio of type-l/type-2 
AGNs at the luminosity of NGC 5548. 

13.4. The Inner Radius of the Torus 

BLR absorption can also help another problem. Bar- 
vainis (1992) was able to satisfactorily explain the <~ 400 
d IR lag in Fairall 9 with dust sublimation temperatures 
of 1300 - 2000 K. However, Fairall 9 lies at the top of the 
IR lag-luminosity relationship shown in Fig. 1 and there 
is an order of magnitude scatter in the lag. NGC 5548 
seems to lie close to the line in Fig. 1. Scaling the IR 
lag of Fairall 9 to NGC 5548 by L 1 / 2 would give a lag 
of ~ 200 d. Since the grains are at ~ 80 Id (or less), 
equation (3) of Barvainis (1987) shows that they would 
be at a temperature of ~ 2500 K. This is too high for 
grain survival. In BLR shielding model, however, the 
continuum seen by the dusty torus in AGNs is absorbed 
by the BLR, and the dust grains can be much closer to 
the central source than they could otherwise be. For a 
factor of 2.7 reduction in the incident energy this allows 
the dust to be at - 120 Id in NGC 5548. ' This is still 
greater than our estimated 80 Id. While this is within 
the observational uncertainties for NGC 5548, there are 
AGNs in Fig. 1 with substantially smaller dust lags for 
their luminosities, so an additional effect is needed. As 
Kishimoto et al. (2007) discuss, dust can survive closer to 
the central source if the dust grains are somewhat than 
those Barvainis (1987, 1992) considered. For larger dust 
grains the IR cmissivity rises substantially (see Draine & 
Lee 1984) so they can cool much more efficiently and have 
lower equilibrium temperatures. There is good evidence 
from X-ray optical observations (Maiolino et al. 2001a, b) 
and optical-UV extinctions curves (Gaskell et al. 2005; 
Gaskell & Benker 2007) for dust grains in AGNs being 
larger than those in the local interstellar medium. 

14. CONCLUSIONS 
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We have used simultaneous observations and the re- 
sults of multi-wavelength variability studies to consider 
in detail the most probable SED of NGC 5548 during the 
highest state in the HST monitoring in 1993. 

We have argued that when the continuum of NGC 5548 
(and other AGNs) is corrected for Galactic and inter- 
nal reddening, detailed photoionization modelling shows 
that the strengths of all the major emission lines can be 
explained by photoionization by the inferred AGN con- 
tinuum so long as the covering factors are large (~ 40%). 
The long-standing "energy-budget problem" can now be 
seen to have been a consequence of not appreciating the 
significant reddening of the continua of most AGNs and 
erroneously assuming a spherical geometry. Because the 
emission-line strengths can be explained by photoion- 
ization from the observed (or extrapolated) continuum, 
there is no need for the additional heating source sug- 
gested by Collin-Souffrin (1986), Joly (1987), and Du- 
mont et al. (1998). 

To reconcile our high covering factor with the lack of 
absorption along the observer's line of sight we need a 
picture of the BLR which is different from the widely 
assumed quasi-spherical distribution of BLR clouds. In- 
stead most of the BLR clouds (especially those far from 
the AGN) have a flattened distribution that extends out 
to the inner edge of the dusty torus. We have shown 
that such a model is consistent with the observed line- 
continuum transfer functions for NGC 5548. 

In the geometry needed to explain the observed line 
strengths, the cloud covering fraction is close to 100% 
in the equatorial plane. We have therefore proposed a 
shielding model of the BLR where the inner clouds ab- 
sorb much of the radiation before it reaches the outer 
clouds or the torus. This solves the problem of LOC 



models not explaining the strong change in cloud condi- 
tions with radius. The shielding BLR picture quantita- 
tively reproduces the ionization-dependent lags seen in 
NGC 5548. 

The absorption by BLR gas also solves several prob- 
lems with the torus. Allowance for BLR absorption in- 
creases the torus covering factor calculated from ener- 
getics, and makes it slightly greater than the mean BLR 
covering factor. This reconciles the torus covering fac- 
tor with the fraction of obscured (type-2) AGNs and ex- 
plains why low-ionization by the BLR is never seen. The 
reduced flux seen by the torus also allows the dust to 
exist at smaller radii and reduces the discrepancy be- 
tween the calculated dust sublimation radii and the radii 
of the inner edge of the torus found from reverberation 
mapping. 

Since the line ratios, equivalent widths, and the con- 
tinuum shape of NGC 5548 are similar to other AGNs 
(see, for example, Sergeev et al. 1999) we expect that 
these results can be generalized to most other AGNs. 
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